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' Q ball solutions are considered within the theory of a complex scalar field with a 

, gauged U(l) symmetry and a parabolic-type potential. In the thin-walled limit, we 

. show explicitly that there is a maximum size for these objects because of the repulsive 

I ■ Coulomb force. The size of Q ball will increase with the decrease of local minimum 

+^ ' of the potential. And when the two minima degenerate, the energy stored within the 

^ I surface of the Q ball becomes significant. Furthermore, we find an analytic expression 

^ ! for gauged Q ball, which is beyond the conventional thin-walled limit. 
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1. INTRODUCTION 

Recently,Theodorakis introduced a piecewise parabolic potential for a complex scalar field (j) and 
he showed that it admits stable Q ball solutions [1]. These solutions have been found analytically, 
unlike the case of Polynomial potentials. The Q ball is a nontopological soliton, and its stability 
depends on whether charge can be lost through emission of charged particles [2,3]. The generation 
of the Q ball[4] and the possibility of phase transition precipitated by soliton synthesis[5] have 
been actively studied. Later work also examined Q-stars [6] and D-stars [7,8] with interesting 
astrophysical implications. The nontopological soliton can appear in a local U(l) invariant theory 
[9-11]. 

Following Coleman [3] , at large distances from the Q ball the field should approach the vacuum 
solution = 0, if the charge is to be finite. It turns out that Q balls exist in a potential 
for large enough values of the charge Q, if the function ^ has a minimum of (p. The simplest 
form of potential is a polynomial one where a negative |^|^ term is needed to make potential dip 
below and a positive |^|^ term to show the potential bounded from below. However, these 

potentials give a difficult task for solving equations of motion. It would be interesting to obtain 
analytic solutions of the Q ball. Fortunately, there is no reason, at least at the classical level, 
to restrict ourselves to polynomial potentials. In fact, Theodorakis [1] has examined a piecewise 
parabolic potential as follows 
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In the 1 < e < 2 case, the global minimum is at = 0, but there is a local minimum at |^| = £, 
the potential being positive everywhere. In this paper, we consider the gauged Q balls where the 
potential is parabalic-type one. Hitherto, ones were restricted to numerical methods or thin-walled 
approximation in the study on gauged Q ball. In the thin-walled approximation, we find that the 
radius of Q ball is increased with the decline of local minimum of the potential. And when the two 
minima are degenerated, the surface term will play an important role in the fuctional of energy. 
Using an iteration method, we have shown an explicit solution of Q ball with gauge interaction 
beyond thin-walled limit in this paper. Furthermore, we can give the asymptotical solutions order 
by order in principle. 



2. QUALITATIVE PROPERTIES OF SOLUTION 

We consider the theory with a complex scalar field 
Lagrangian density is 



coupled to a U(l) gauge field A^^. The 



U 



(2) 



where F^i^ = dfj,A^ — d^A^^ and C/(|0|) is the parabolic-type potential (l).The conserved current 
and the energy momentum tensor which induce the conserved quantities are 
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(3) 



and 



T — - 



(df, - ieA^)(l){dy + ieA^)4> + {dy - ieAy)4>{di^ + ieA^)<j) 
respectively. To find the solutions of Q ball, we use ansatz as follows 



(4) 



iuit 



and 



Af, = Ao{r)Sf,o 



(5) 



(6) 



where we assume a; > for definiteness. The lowest energy state will have no electric currents and 
therefore no magnetic fields. The spatial components of gauge potential are zero as there is no 
magnetic field. We choose a boundary condition Aq ^ as r ^ oo. 
Using the above configurations, the Lagrangian becomes 



L = Att drr^ 
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(7) 



where g = lo — eAQ{r), and prime denotes d/dr. The Noether charge associated with the U(l) 
symmetry becomes 



Q = j d^xgf 



(8) 



When we insert Eqs. (5) and (6) in the general field equation deduced from Eq. (2) we find the 
following equations of motion: 



f" + -f' + fg' = f-e, 

r 



f" + -f' + f9' = f, 
r 



9" + -g' = e'fg . 
r 



r<R,f>l 



r>R,f<l 



(9) 



(10) 



(11) 



The fact that the field equations (9-11) are nonlinear implies that the analytic solutions that will 
be found will be approximate. To gain insight into the Q ball for e 7^ 0, we show some qualitative 
properties of the solution. Prom (8) and (11) we have 



g{r) u 



(12) 
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as r — > oo. Then Eq. (10) can be reduced to the asymptotic form 

l^rfr+iu;^-l)f = (13) 
r 

For the approximate equation (13), we have a solution as fohows 

/(r) = foex.p{-Vl-u;'^r)/r (14) 

Clearly, a necessary condition for the existence of a solution is a; < 1. Additionally, for the solution 
to be well-behaved at the origin, /' and g' should approach zero at least faster than r for r ^ 0. 
By using the boundary conditions of / and g the energy functional can be written as 

E = ^uQ + A7rJ r^dri^f'^ + V{f)] (15) 

For the e 7^ case, we expect that the energy will be increased over the e = case due to Coulomb 
repulsion, with Coulomb energy becoming more important as Q becomes large. Therefore, if 
dE/dQ > 1, we must consider that some charge can be put into the interior region of the Q ball 
and some can be put in free particles. As discussed above, there exists a Qmin and a Qmax, so that 
when Q > Qmin the soliton is quantum mechanically stable and when Q > Qmax the lowest energy 
state of the system is composed of a soliton with charge Qmax together with free particles carrying 
charge Q — Qmax- In the e = theory, since Qmax 00, the condition Qmin < Qmax is always 
satisfied. 



3. THIN- WALLED APPROXIMATION 

If we choose a continuous ansdtz for the scalar field, the most convenient one might be the 
piecewise ansdtz 



fir) 



f r<R 
f-i{r-R) R<r<R + S (16) 
r > R + S 



where 6 is the width of the wall and / is a constant. But unfortunately, this ansate,unlike e = 
case, does not satisfy the equations of motion of the scalar field. Therefore, one must choose 
the thin-walled approximation in which the spatial configuration of the scalar field /(r) can be 
considered as a step function. 

In the thin- walled approximation, the spatial configuration of the scalar field /(r) can be con- 
sidered as a step function 

f{r)=~fe{R-r) (17) 

Although the gradient of / in the thin surface is a delta function, and hence the associated density 
of kinetic energy is infinite, the surface term in the integral of energy is finite. Especially,in the 
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case of e < 2, the contribution of the kinetic energy term to the integral of energy is neghgible. 
Therefore the constant / in Eq. (17) is determined by minimizing energy functional. In this case, 
we can solve Eq. (11) 



(.;-^)i?sinh(r/r) 
sinh(e/i?) r 



9{r) 



g{r) = u 



for r < R 
for r > R 



(18) 
(19) 



where the gauge is chosen so that — for r — oo. Inserting the above solution into Eq.(8), we 
have 



UJ = 



AttR 



1 



tanh(e/i?) 
efR 



The total energy of the ball is reduced to 

e'^Q^l^ tanh(e/i?)^-^ 



SttR 



1 



efR 



+ fuif)R^ 



(20) 



(21) 



Prom dE/dR 



= 0, the radius R* of the ball is determined by 



^T,^Jw{f)/T^R 



efR 



1 



e^Q 



(22) 



.tanh(e/i?) 

It is a formidable task to find exact solutions of the transcendental equation (22) . We have to use 
e/i? <C 1 approximation [9]. Eq.(21) can be reduced to the following form 



E = Q 



\'^u{f)^ 


± 

2 


y p \ 





1 + 



(|f/Q)^(/V2^(/))^ 



(23) 



Minimizing the energy with respect to / and Q{R) we find that 

/* = 2- 
3Q 



5e 



R* 



Btt V2(2 - e) 



1 + 



Stt 

1 



3e^Q 



and 



E^ = Q\ 1 



45 V7rv'2(2-£ 



(24) 
(25) 



£ r 1 

21 5 V7rv/2(2-£) 



(26) 



For a fixed charge Q, we see that both the radius and the energy are larger than ungauged one. 



Furthermore, from the condition dE/dQ 
charge Qmax as 



Q=Qr, 



1, one can obtain the upper bound on the total 
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Q - — 



(27) 



When Q > Qmax the lowest energy state of the system is composed of a Q ball with charge Qmax 
together with free particles carrying charge Q — Qmax- Here note that the inconsistency of the 
approximation efR <^ 1 when e approaches 2. In the £ = 2 case,we must considered the surface 
term with a finite width. If efR S> 1, as will be discussed later, the energy can be approximately 
written by 



e2n2 

E = + AireR^ (28) 



From dE/dR 



= 0, the radius of the Q ball is given by 



and 



The upper bound on charge Qmax 

is also given by 



E^ = E{R^) — ^ ^\ 



(30) 

237r3 



TT 

(31) 



For simplicity, we write /*, i?* and E^ as f , R and E hereafter. 



4. BEYOND THIN- WALLED LIMIT 

We may use the successive approximation method of differential equation for higher order 
approximation of the Q ball solution. Using rg(r) = e"^''^ and v(r) = u'(r), Eq.(ll) can be reduced 
to 

v' = e^f - (32) 
with the initial condition v = vq for r = ro- It can be written in the form 

v = vo+ r{e'^f-v'^)dr (33) 

Jro 

Substituting the thin- walled solution fi(r) instead of v to the right member, wc obtain a new 
function V2 different vi, unless vi is a exact solution of Eq.(32). Substituting V2 instead of v to 
the right member of Eq. (32), we obtain a function 113 and so on. The sequence vi, V2, ■ ■ •, Vn, 
■ ■ ■ obtained in this way is convergent in a certain interval containing tq to the desired solution of 
the given equation, provided that the assumptions of Cauchy's existence theorem. In principle, we 
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can find the exact solution g[r) if /(r) has been known. This method is also called the iteration 
method. 

In the weak coupling situation e <C 1 and g{r) = uj, we have solutions of Eq. (9) and (10) as 
follows 

^^^^^sinhM (34) 
r 

/ = ^e-'^R-r) r>R (35) 

where v = \/\ — uS^^ = 17^ ^'^d 6 = (1 — ^) sinhz/R - -^"-"^ thin- walled limit, we have / = a + in 
the Eq. (19). Equivalently, we have 

-i;i(r) = e(a + 6z/)coth[e(a + 6z/)r] r<R (36) 

vx{r) = (r r>R (37) 

Substituting v\(r) instead of v to Eq. (33), we obtain high-order approximation 

-|-e(a 6z/){coth[e(a + 6i/)r] - e(a + 6zy)r} + C r <R 



^;2 = [r---^\ -1^1^ 2e^R'^i'eyi^{2vR)Ei{-2vr) r>R (38) 

V Attu! J r 

where Shi{x) is the hyperbolic-sine- integral function and Ei{x) the exponential-integral function. 
Using the relation 

Ei{-2ur) = -e^'^' % (39) 

^ ' Ji 2ur + lntt^ ^ ' 

one can easily find that v{r) satisfies the boundary condition at space infinity v{oo) = 0. The 
integral constant C can be fixed by continuity at r = -R. Therefore, we have 

/ e'^Q\ sinh\e(a + buh-] / r\~ ^ .e^ , , .9 9w„9 9n 
52 = w---^ — TI-7 e^^{ — {2ahv + h^v^){R^-r^) 

2 

+2e^ah[rShi{vr) - RShi{vR)\ + -e^ah[cos\i{vR) - cosh(i^r)] 
+e%'^v[rShi{2vr) - RShi{2vR)] + e%'^[cosh{2vR) - cosh(2i/r)] 

WR^v + 2e^R^veyLY>{2vR) Ei{-2vr)dr} r<R (40) 

Jr 



/oo 
Ei{-2ur)dr]} r> R (41) 
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Eq. (34), (35), (40) and (41) are the solution of high-order approximation beyond the thin-walled 
limit. Our approximate analytic solution is well in agreement with the numerical results [9], which 
have been obtained by many authors. The study on Q ball in a parabolic potential, as a toy 
model, can enable us to understand the Q ball analytically, which provide more information than 
numerical study. In principle, by continuing the procedures mentioned above, we can obtain the 
approximate solution with any high-order accuracy. It's obvious that this method can also be 
employed in investigating other spherically symmetric soliton, and the corresponding results will 
appear elsewhere. 
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